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Abstract. The penetration, diffusion and slowing down of electrons in a semi- 
infinite air medium has been studied by the Monte Carlo method. The results 
are applicable to the atmosphere at altitudes up' to ~ JOO^km. Most of the 
results pertain to monoenerge tic electron beams injected into the atmosphere 
at a height of 300 km, ei:her vertically downwards or with a pitch-angle distri- 
bution isotropic over the downward hemisphere. Some results were also obtained for 
various initial pitch angles between 0° and 90°. Information has been generated 
concerning the following topics: (a) the backscattering of electrons from the 

atmosphere, expressed in terms of backscattering coefficients, angular distri- 
butions and energy spectra of reflected electrons, for incident energies T Q 
between 2 keV and 2 MeV; (b) energy deposition by electrons as a function of the 
altitude, down to - 80 km, for T^ between 2 keV and 2 MeV; (c) the corresponding 
energy, deposition by alectron-produce‘3 brcmsstrr^lung , .flown***" ~ ™ Vm; (d) rim 
evolution of the electron flux spectrum as function of the atmospheric depth, 
for T between 2 keV and 20 keV. As far as possible, the results have been 
expressed in a scaled form which reduces the explicit dependence on, and permits 
interpolation with respect to, T Q . Energy deposition results are given for 
incident electron beams with exponential and power-exponential spectra. 
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1 . INTRODUCTION 


In the course of earlier Monte Carlo calculations of auroral luminosity 
patterns (BERGER, SELTZER and MAEDA, 1970; to be referred to as BSM) and breros- 
strahlung flux spectra in the -atmosphere (BERGER and SELTZER, 1972) we have 
obtained other electron and bremss trahlung transport results which may be useful 
for the analysis and interpretation of observed atmospheric phenomena. These 
results are given here, and pertain to the following topics: (a) the backscat- 

.tering of electrons from the atmosphere, as 'the result of multiple Coulomb 
scattering; (b) the altitude-dependence of energy deposition by electrons pre- 
cipitated into the atmosphere; (c) the corresponding altitude- dependence of energy 
deposition by electron-produced bremss trahlung ; and (d) the modification of the 
energy spectrum of electron beams as they. penetrate down into the atmosphere. 

2. MONTE CARLO MODELS 

The Monte Car 1 ^ method involve z the simulation of -all important physical 
processes by random sampling. One must sample the elastic scattering of electrons 
by atoms, the inelastic scattering of electrons by. orbital electrons, the pro- 
duction of bremss trahlung, and the Compton scattering and photoelectric absorption 
of the bremss trahlung photons. 

The number of elastic scatterings which an electron undergoes before it is 
slowed down to rest in the atmosphere is generally . quite large. For example, with 
the cross sections used in BSM, one finds that the average number of collisions is 170 
for an initial energy of io keV , 290 for 20 keV and 1060 for 100 keV. An alterna- 
tive, less time-consuming, method has therefore been used in most of the calculations 
reported here, which is referred to as Monte Carlo Model A in BSM. This approach 
reduces the required computational effort through the combined use of random sampling 
and analytical multiple scattering theories. Each electron track to be sampled is 
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divided into a number of segments (usually not more than ~ 100) whose lengths are 
chosen so that the number of elastic scatterings per segment is large whereas the net 
multiple-scattering deflection and energy loss per segment are small.*/ The angular 


, In the present work the lengths of these segments were chosen so that - on the 
average - the energy loss per segment was equal to 4.2% of the electron's energy at 
the beginning of the segment. The angular multiple-scattering deflection was 
sampled at the end of each segment. The actual energy loss was sampled at the end 
of every second segment. Energy-dependent cross sections at intermediate points 
along the trajectory were evaluated assuming a linear dependence of energy on path- 
length traveled. 

deflections in successive segments are sampled from the distribution of GOUDSMIT and 
SAUNDEPvSON (1940), and the energy-losses from the distribution of LANDAU (1944) with 
the binding correction of BLUNCK and LEISEGANG (1950),. These deflections and energy 
losses are then combined to construct the complete electron track. Monte Carlo Model 
A also includes the production and subsequent diffusion of secondary electrons, and 
the production, scattering and photo-electric absorption of secondary bremsstrahlung 
photons. The transport of bremsstrahlung photons is followed by conventional random 
sampling techniques (FANO, SPENCER and BERGER, 1959) involving the sampling of all 
successive Compton scatterings until photoelectric absorption occurs. 

In order to obtain our backscattering results for initial electron energies 
^ 20 keV we have used another method, which is described in detail in BSM where it 
is referred to as Monte Carlo Model B. In this model all successive elastic scat- 
terings of an electron are followed by random sampling, with use of the appropriate 
single-scattering cross section. The numerous inelastic collisions are lumped 
together and are treated in the continuous-slowing-down approximation, assuming a 
relation between electron energy and pathlength traveled given by stopping power 
theory. ■ In this approximation, secondary knock-on electrons are not taken into 
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account. 


*/ 


In a radiobiological calculation for a water medium (BERGER, 1972) a more elaborate 
method has been introduced which could be designated as Model C, In this model all 
hard inelastic collisions are sampled individually, whereas the soft inelastic 
collisions are treated in the continuous- si owing- down approximation with use of a 
’'restricted stopping power" theory. Hard inelastic collisions are defined to be 
those in which knock-on electrons with an energy greater than a 'chosen threshold 
value, say 200 eV, are set in motion. This procedure makes allowance for energy- 
loss straggling and the subsequent transport of the knock-on electrons. 


For. .each source energy and geometry of interest, a large set of electron tracks 
(typically 10,000) was simulated according to Model A or B or, in some cases, with 
both. The sampled tracks were then analyzed, to obtain information about various 

transport phenomena including backseat tering , energy deposition and electron flux 
spectra. Details of some of the sampling procedures and analysis are described in 
BERGER (1963) and BERGER and SELTZER (1968). 

3. CROSS SECTIONS 

In Monte Carlo Model A the cross section for the elastic scattering of electrons 
by atoms, needed for the evaluation of the Go udsra-it-. Saunders on multiple-scattering 
distribution, is taken from the theory of MOTT (1929) which includes spin and 
relativistic effects. In Model B, applied at low energies, relativistic effects 
are unimportant and the Rutherford cross section can be used. In both models, the 
cross section at small angles is modified by a. screening correction given by 
MOLIE^E (1947, 1948). The electron stopping power values have been taken from the 
Bethe theory as formulated by ROHRLICH and CARLSON (1954). 

Because of the use of the Born approximation, and because of the incomplete 
treatment of atomic binding effects the .theore tical cross sections used in our 
Monte Carlo program to describe interactions between electrons and the atmospheric 
constituents are expected to be applicable only down to some cut-off energy, perhaps 
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comparable to the K-shell binding energies of the atoms involved. estimate that 
the overall cross section uncertainty is still only 10% to 15% at 0.5 keV, but 
becomes greater at lower energies. We expect that the transport calculations of 
energy deposition will give reliable results under conditions such that the electrons 
have lost the predominant part of their initial energy t and have traveled a path- 
length equal to a large fraction of their initial range, when they reach an energy 
of 0.5 keV. . 

- There is experimental evidence regarding the reliability of the calculations. 

It has been shown in BSM that the calculated spatial distribution of energy 
deposition in air is in good agreement with the results of a laboratory experiment 
by GRtJN (1957) with electron beams incident with energies as low as 5 keV. Com- 
parisons with similar measurements of COHN and CALEDONIA (1970) in nitrogen indi- 
cate similarly good agreement down to 2 keV. In this paper it will be shown that 
the calculated backsca ttering coefficients for air are in agreement with experi- 
mental data down to 2 keV. Thor** are at present no laboratory data which could be 
used to check the calculated electron flux spectra. It is possible that the spectral 
distributions are more sensitive to the cross section input than energy deposition 
distributions or backsca ttering coefficients. The electron flux spectra are there- 
fore expected to be the most tentative of the results in this paper, and subject to 

* / 

revision in the light of better cross section information. ' 


New information on low-energy electron scattering and energy loss processes is 
gradually accumulating. We are beginning a new cycle of transport calculations 
according to Monte Carlo Model C, using data such as the semi-empirical energy 
loss functions of GREEN and PETERSON (1968), the inelastic electron scattering 
cross sections measured by OPAL, BEATY' and PETERSON (1972) and calculated by 
OMIDVAR, KYLE and SULLIVAN (1970), and the measured elastic scattering cross 
sections of BROMBERG (1970) and SHYN, STOLARSKI and CARIGNAN (1972). Much work 
is yet needed to generate from these and similar data the comprehensive data base 
needed for the Monte Carlo calculations. 
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4. SCHEMATIZATIQN 


Electrons are assumed to be incident onto the boundary of a semi-infinite 
air medium with kinetic energy T q and obliquity 0 q (angle between the direction of 
incidence and the normal to the boundary plane). Electron transport is calcu- . 

lated as a function of the depth z in the medium, i.e, the mass thickness (in 

-2 

g cm ) measured from the boundary plane. The boundary plane is treated as non- 
electron— pto duce d 

reentrant, so that electrons (or / bremsstrahlung photons), after escaping 

from the medium, cannot return to it. The incident electron beam current across 

the boundary plane is assumed to be either monodirectional (with a fixed obliquity 

angle 0 Q ), or characterized by an angular distribution proportional to cosQ^. 

In applications to the atmosphere, the obliquity angle 0^ (at sufficiently 

high magnetic latitudes) can be identified with the incident pitch angle; z is 

m 

pitch-angle 

the atmospheric depth; the cosine-law/distribution corresponds to .the case of an 
incident electron flux isotropic over the downward hemisphere (to be denoted by the 
abbreviation IDH) . As has been shown by REES (1964a), the relation between the 
atmospheric depth z^ and the altitude h is given by 

« 1/2 

= J D < h > [frill* ■ ’ a) 


m 


m 


where 


x = cos X (1 + h/a) , 
m ’ 


( 2 ) 


is the geomagnetic latitude, a is the radius of ’the earth and D(h) is the 
density of the atmosphere. In the present work we have assumed that D(h) is given 
by the CIRA (1965) mean atmosphere for heights between 30 km and 300 km. 


5. NEGLECT OF MAGNETIC MIRRORING 


Electrons are deflected in the atmosphere not only as the result of multiple 
Coulomb scattering but also by magnetic mirroring. The combined action of both 
effects' has been included in the calculations of WALT, Me DONALD and FRANCIS (1967), 
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who solved the electron transport equation in the Fokker-Planck approximation, and 
in the Monte Carlo calculations of WEDDE (1970) and McENTIRE (1972). The calcu- 
latioris reported here do not 'include the effects of magnetic mirroring. The re- 
sults are nevertheless applicable to atmospheric problems, provided one limits the 
applications to altitude no greater than, say 300 Km. * have therefore assumed 
that the atmospheric depth ^ = 0 (the boundary of the semi-infinite air medium 
where the characteristics of the incident electron beam are specified) corresponds 

to 311 altitude of 300 1cm. 

This assumption has practically no effect on the Interpretation of the calcu- 
lated multiple Coulomb scattering results, because the mass thickness of the atmosphere 
above 300 km is negligibly small compared to the range of the electrons at the 
energies of interest. The magnetic deflection angles ire generally small compared 
to Coulomb multiple-scattering deflections for electrons starting out at 300 km. 

They can be estimated assuming a magnetic dipole field and using the adiabatic 
invariance of the ^entity sin 2 8/H, where 6 is *he P- Uh angle and H the magnetic 
field strength. Epical pitch-angle changes are given in Table 1 for electrons 
that start with various initial pitch angles at altitudes of 300 or 200 km and 
travel down to heights of 90 or 50 km. It can be seen that in most cases the 
magnetic deflection angles amount to only a few degrees and become significantly 
greater only when the initial pitch angle approaches 70°. Moreover, the fraction 
of the electrons with pitch angles greater than 70° is small for an incident 
electron flux Isotropic over the downward hemisphere. 


6. ELECTRON BACKSCATTERING 

The backscattering of electrons from the ionosphere was observed in rocket 
experiments by -McDIARMID, ROSE and BUDZINSKI (1961) and in the Injun III satellite 
experiments of O'BRIEN (1964). According to O'Brien, approximately 107 o of the 
electrons reaching the atmosphere are backscattere d into the magnetosphere, which 
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is consistent with the results of our calculations. More evidence on backscattering 
has been obtained in an "electron echo" experiment by HENDRIKSON, McENTIRE and 
WINKLER (1970) in which 40-keV electrons were injected into the atmosphere from 
a rocket-borne accelerator at. a height of 300 km at a geomagnetic latitude of 
52° N, with emission pitch angles between 60° and 120°. The injected electrons 
traveled back and forth between the mirror point near the rocket and the conjugate 
southern mirror point, and several electron echos were detected for each emitted 
pul-se. The southern conjugate mirror point was almost at sea level, so that 
Coulomb scattering rather than magnetic mirroring was primarily responsible for the 
appearance of the echos. In an analysis of this experiment, McENTIRE (1972) has 
made an elaborate Monte Carlo calculation which takes into account not only multiple 
Coulomb scattering in the atmosphere (by a method similar to our Monte Carlo Model B) 
but also the motion of the electrons under the action of the geomagnetic field. 

These calculations are reported to account quite well for many features of the 
observed electron ^ch^q. 

Our calculations are aimed not so much toward the explanation of a particular 

experiment as toward the production of background information of general interest. 

First of all, we shall examine the validity of the multiple- sea ttering calculations 

through the comparison with laboratory experiments . Such a comparison is possible 

in respect to the number albedo R^, i.e. the fraction of the incident electrons that 

are backseat tere d from a semi-infinite medium.' Several authors have measured the 

number albedo for series of target materials with atomic numbers bracketing that of 

air (z = 7.4). We have interpolated graphically with respect to atomic number and 

have thus obtained the experimental R^- values shown in Fig. 1. Also shown in that 

figure is the calculated curve of R vs. the energy of incidence. T for T between 

N o o 

2 keV and 2 MeV. The calculated results have an estimated statistical uncertainty 
(relative standard deviation) of 37.; including systematic error, the uncertainty 
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may be as large as 5%. The uncertainty of the experimental results is indicated 
in some cases, and can also be judged from the dispersion among the backscatter ing 
coefficients from different experiments. Within the combined limits of experi- 
mental and calcula t iona 1 error there seems to be good agreement between predicted 
and observed number-albedo values. The theoretical number albedo shows a systema- 
tic increase as is decreased. So do the experimental values, except that some 
experiments indicate a leveling off at an energy T 2 keV.* 

•ft 

Calculated backscatter ing coefficients for air have been given previously by 
MAEDA (1965). His results, for T q = 20 keV, are in good agreement with the preser.t 
results. At lower energies his backsca ttering coefficients are smaller and have a 
different energy dependence, showing a decrease rather than an increase with de- 
creasing T These discrepancies may be related to difficulties with Maeda's 
electron flux calculations discussed below in Section. 9. 

The calculation also provided results not available from experiments, namely, 
the number albedo and the energy albedo (fraction of incident energy back- 
scaLi.ci.ed) as functions of the obliquity of the incident electrons. In Fig. 1,1^ 
is shown together with for incident energies between 2 keV and 2000 keV, for 
perpendicular incidence (9 q = 0 ) and for a cosine-law source (IDH case in the 
atmosphere). Figs. 2a and 2b show and R £ as functions of the incident obliquity 
angle 9 q , for incident energies of 10, 100 and 1000 keV. 

have insufficient data to give a distribution of backsca ttere d electrons 
differential in energy and direction, ye can only show the angular distribution 
regardless of energy, and the energy spectrum regardless of the direction of 
emergence. The dependence of these single-variable distributions on the incident 
energy T q is greatly reduced if one normalizes them through division by the back- 
scattering coefficient Normalized angular distributions w(9) are show in 

Fig. 3 for G o = 0 and for the IDH case, in the form of universal curves which are 
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applicable for all T -values between 2 keV and 2 MeV. .Normalized cumulative 
energy spectra are shown in Figs. 4a and 4b for 9^ = 0 and the IDH case. Again 
the results are insensitive to the value of T , particularly for < 100 keV, 

The dependence of the normalized cumulative energy spectra of backsca ttered 
electrons on the incident obliquity angle is shewn in Figs. 5a, b and c for 
incident energies of 10, 100 and 1000 keV. 

7. ’ENERGY DEPOSITION’ BY ELECTRONS 

The basic quantity of interest is the energy depos it ion function A(z^) which 

is defined as the energy deposited per unit mass thickness at a depth from 

the boundary of a semi-infinite medium. We shall assume that A(z^) is expressed 

in units of eV/(g cm -2 ), and the depth in g cm -2 . The amount of energy deposited 

per cm 3 in the atmosphere at an altitude h, normalized to an incident beam of 1 

electron cm 2 , is given by D(h) ACh(z^) 1 , where D(h) ?.s the atmospheric density, 

in e cm -3 and where z and h are related by Eq(l). We shall refer to the product 
° ’ m 

D(h)A[h(z^)] as the altitude profile of energy deposition . 

*Th e notation is the same as in BSM where the symbol p was used to indicate a 
radial spatial variable and D to indicate the density. The energy deposition 
function A was called penetration function in BSM. 

The energy deposition function A(z^) has' previously been computed in BSM 
with Monte Carlo Model B for electrons with incident energies between 2 keV and 
20 keV. Using Monte Carlo Model A, we have now extended these calculations up 
to 2000 keV for the case of perpendicular incidence, and to 500 keV for the 
cosine-law (IDH) case. At 20 keV, the old and the new calculations are in close 
agreement. The accuracy of the calculated energy deposition functions is esti- 
mated to be 37 0 or better. As shown In BSM, there is good agreement between the 
calculations and experimental laboratory data in air obtained by GRlJN (1957) at . . 
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energies up to 54 keV. Our results are therefore also consistent with the empiri 
cal energy deposition/of n REES n (1963) at energies up to 300 keV, obtained by interpo- 
lation and extrapolation of Grun's data. Further experimental confirmation is 

available from the measurements of McLAUGHLIN and HUSSMANN (1969) at 100 and 400 
keV, and the measurements of R0SENSTEIN, EISEN and SILVERMAN (1971) at 2000 keV, 
in a medium (polystyrene) with an average atomic number close to that of. air. 

In order to minimise the dependence of the energy deposition function on the 
initial electron energy T 0 , it is useful to tabulate it in scaled form, by ex- 
pressing depths as fractions of the initial electron range. The simplest scaling 
parameter that can be used is the mean range* 7 r Q , calculated by integrating 


*The mean range r (sometimes also denoted as c.s.d.a. 
lated in the continuous-slowing, down approximation) is 
pathlength which the electron would travel if it lost 
its trajectory at a rate equal to the mean loss given 
(see, e.g. BERGER and SELTZER, 1964). 


range because it is calcu- 
numerically equal to the 
energy at every point along 
by stopping-power theory 


the reciprocal of the stopping power from energy T q down to zero. For low values 

of T say below 5 keV, the mean range r Q is not very well defined because of lack 
o * 

of adequate information about low-energy cross sections. In BSM this difficulty 

a **/ 

was circumvented by using as/scaling parameter the so-called practical range r p 


**The practical range r is the quantity obtained by extrapolating the linear 
portion of the curve Bf energy deposition vs., depth (see BSM) ; r depends no 
only on the initial electron energy, but also on the distribution of 
obliquity angles, and is thus different for perpendicular incidence and the IDH 

case . 


which is better defined at low energies but must be obtained from a complete trans- 
port calculation. In the present work, for T q 2 20 keV, we have reverted to the use 

of the mean range r . 

The scaled energy deposition function will be found tabulated in Table 2a for 
the case of perpendicular . incidence, and in Table 2b for the IDH case. Each table 
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contains a high-energy part, consisting of. the dimensionless quantity ( r 0 /T Q ) A(z^) 

as a function of z /r , for T S 20 keV; and a low-enfergy part, excerpted from BSM, 

mo o 

consisting of the dimensionless quantity (r /T ) A(z ) as a function of z /r , 

p o m m p 

for T q £ 20 keV. Together, the low and high-energy parts span the energy region. 

from 2 keV to 2000 keV for perpendicular incidence, and from 2 keV to 500 keV for 

the IDH case. At the matching energy (T q = 20 keV) the two different scaled 

representations give the same values for A(z ). The values of the ranges r and • 

m o 

r^, needed to interpret the data in Tables 2a and 2b, are given in Table 3. 
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.. t-o T one can readily obtain energy deposition 
By interpolation with respec 0 

f , of energy deposition for any incident electron 

functions and altitude profiles of energy cep 

and power-exponential spectra of 
spectrum of interest, for example, exponent P 

r 2b- and 3 are given only for incident 

the form T f exp(-T /«>. The data in Tables 2a, 

o £ interest may contain 

energies T greater than 2 keV. However, the spectr 
" 8 ° . ^ nif leant contribution to the energy 

electrons with lower energies that make a sign 

deposition, particularly at very high altitudes. It is possible to inc « 
contribution from these electrons by extrapolating the scaled energy epos^ 
function, (r /*,) *(..), the practical range. r p . to lower energies. 

t J, : to assume that the scaled energy deposition function retains 

ri: - . «. — ~ - — — - - 

u . tlon is probably smaller than 20 X even for T -valu-s 

error incurred by this assumption is prooao y 

, . - _ for the IDE case leads to estimated values 

as low as 0.2 keV. Extrapolation o£ i p f . 

-6 . n -6 . , 6 x 10 - 7 g cm' 2 at 1, 0.5 and 0.2 keV, respectively, 

of 4.7 x 10 6 , 1.8 x 10 . and 6.6 x 10 g . 

. . . -„ r ve have assumed .that -the ex r -nen ia 

When calculating tne given -el 

or power-exponential' spectra extend down to T q = 0. The contributio 

' trons with initial energies less than 2 keV to.the energy deposition are shown in 

, -pie 6 ve give the altitude profile 

Table 4 for various atmospheric depths. In Fig. S 

of energy deposition for a purely exponential spectrum 0)» assuming an 

incident flux isotropic over the downward hemisphere and e-folding energies a 

between 5 keV and 200 keV. We estimate that the uncertainty of these results 

a - • <- 300 km £ 10% at 200 km f and ~ 

due to the low-energy extrapolation is £ 20/. at 300 km, ~ 

at 150 km. 
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In Table 5 we give data for the case of a power-exponential spectrum 
exp(-T/a): the height h at which the altitude profile peaks and the cor- 
responding peak value D(h)A[z (h)], for various a-y combinations. All of the 

m 

above results are evaluated assuming in Eq.(l) a geomagnetic latitude of 90°. 

The dependence of the altitude profile on the geomagnetic latitude is slight, 
and is indicated in Table 6 for the example of a power-exponential spectrum 
with y = 0.5 and a = 10 keV. 

8. ENERGY DEPOSITION BY ELECTRON-PRODUCED BREMSSTRAHLUNG 

Electron beams deposit their energy in the atmosphere not only directly 
but also indirectly via bremss trahlung . The amount of bremss trahlung energy 
deposited is small but significant because the p ho tens can pene- 
trate down to altitudes of 30 km- or lower, which the primarily electrons cannot 
reach. 

•.The energy -d eposition - funculon^ f or- >. l^em& sx' rah Iung , -(z ) has been obtained, 

by evaluating the expression 

T o 

A (z ) e [ 4 (k,z , T ) kp (k)dk (3) 

or mJo m oen 

o 

where k is the bremss trahlung photon energy, 0 (k,z ,T ) the bremss trahlung flux 

■q m o 

spectrum at depth z due to incident electrons of energy T , and u (k) the 

m o en 

photon energy absorption coefficient for air (HUBBELL, 1969). The bremss trahlung 
flux spectra were taken from recent calculations (BERGER and SELTZER, 1972) for 


13 



the case of broad beams of electrons entering a semi-infinite medium with a 
cosine-law angular distribution (corresponding to the IDH case and wide-area 
precipitation into the atmosphere). In order to reduce the explicit dependence 
0n V we k ave it useful to tabulate the quantity as a function 

of the variable z^T^. Such data are given in Table 7 and are sufficient to 

determine, by interpolation, the energy deposition by brems strahlung at atmospheric 

-2 

depths up to 20 g cm , for incident electron energies between 20 keV and 2000 keV. 
By such an interpolation we have obtained altitude profiles of energy deposition 
by brems strahlung resulting from the uniform wide-area precipitation into the 
atmosphere of electron beams with exponential energy spectra. These results, for 
e-folding energies between 5 keV and 200 keV, are shown in Fig, 7. 


9. ELECTRON FLUX SPECTRA. 


. As an electron beam penetrates down into the atmosphere, the spectrum of 
primary electrons.. Is. ..shifted .towards .lower energies, -and a -buildup of lo w-enp-rgy 
secondary electrons takes place. Simple estimates of the change of the primary 
spectrum have been made by various authors, e.g. -REES ( 1964b ) and ST0LARSKI (1968) 
who implicitly assumed that the spectrum at any depth z^ is monoenergetic and con- 
centrated at an energy 


m 


T ^ Z m^ T o ' J A ^ z, m) dz ' m > 


(b) 


where A(z,^) is w he energy deposition function. This is equivalent to assuming 
that the electrons move along straight lines and lose energy at a rate per unit 
pathlength equal to A(z ffl ). 

The spectrum actually is not monoenergetic but rather broad 

because of the occurrence of two types of fluctuations. 
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First, the electrons do not really travel in straight lines, because of multiple- 
scattering angular deflections, so that the pathlength traveled down to depth z ffi 
is a stochastic quantity. 'Second, even for a given pathlength, the actual 
energy loss fluctuates around the mean given by the continuous -slowing -down 
approximation. Both types of fluctuations - pathlength and energy-loss straggling 
have been taken into account in the present calculations which also include the 
buildup of secondary electrons. 

* We have calculated the spectrum only down to a cut-off energy A, which was 
chosen to be % of the initial energy T q , and have also computed the average. 

* . . 

For the lowest source energy treated, 2 keV, the cut-off energy A was 0.1 T = 

200 eV. * ~ ° 


number of electrons per unit depth, n‘ A (z ), that reach an energy A . The relation 
between the flu:-: spectrum F(T, z ) and the energy deposition function A(z m ) is 
given by 


A (z ) = 

m 




<o 


nr 


F(T,Z ffl )L A (T)dT 


( 5 ) 


where F(T, z ), in MaV is the spectrum of primary plus secondary electrons, and 


The electron flux differential in energy and direction is defined as the number 
of electrons per unit energy and solid angle that cross a unit area perpendicular 
to the direction of momion. The spectrum F(T, z ) is obtained from the doubly dif- 
ferential flux distribution by integrating over all directions of motion, and can 
thus be interpreted as the distribution with: respect to energy of electrons cross- 
ing a small spherical probe at depth z .. 


where L a (T) is the restricted stopping power. L^(T) includes energy losses from 
electron interactions resulting in atomic excitations or in ionizing events in 
which secondary electrons with energy smaller than A are released; it is less than 
or equal to the ordinary stopping power and can be computed with the use of 
Eqs(22-23) of BERGER and SELTZER ( 196 A). 
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Comparison with a Previous Calculation , A Monte Carlo calculation of electron 
flux spectra, somewhat similar in scope to the present work, was done by MAEDA 
(1965) whose results were fitted to analytical expressions by MAEDA and AIKIN 
(1968) and used by REES (1969) for the analysis of satellite and rocket measure- 
ments of auroral electrons. Recently, SHEMANSKY, DONAHUE and ZIPF ( 1972 ) have 
studied the relation between the emission lines from excited nitrogen molecules 
and the energy flux spectrum of the auroral electron responsible for the excita- 
tion. Their analysis led them to the conclusion that Maeda's spectra are deficient 

in low-energy degraded primary electrons, and that the rate at which electrons 

five 

from these spectra would ionize the atmosphere is up to / times smaller than the 

integrated ionization rate inherent in the incident electron flux. The present 

result supports the conclusions of Shemansky et al . 
a 

For/ detailed comparison we have chosen the case of 20-keV electrons incident 

.perpendicularly. ..Fig, S shows .the .electron spectrum R(X,z^) calculated with our 

Monte Carlo Model A and the corresponding quantity, denoted as "dif ferential energy 

spectrum'* given earlier by MAEDA (1965). In Fig.. 9 we show the distribution function 

n A (z ) also obtained in our Monte Carlo calculation. "When F(T,z ) and n.(z ) are 
Am m A m 

inserted into the right-hand side of Eq(5), the energy deposition thus obtained is 

in agreement to within 5% or better with the value of the energy deposition function 

•k f 

A(z^) from Table 2a. This shows that our energy deposition and flux spectrum 
•k 

In the case under discussion, the contribution to energy deposition from electrons 
with energies £ A (not treated explicitly in the Monte Carlo calculation) is 

relatively minor, and can be approximated by the formula 

[A-n.(z )/A(z )1 = 0.064 + 0.122 (z /r ) . 

A m m m o 

calculations are indeed consistent in the sense that energy is conserved. 
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It can be seen from Fig. 8 that Mae da* s "differential spectrum" is lower than 

F(T z ) at all three depths shown, and that the discrepancy becomes more marked 

m */ 
at the lower spectral energies. When a weighted integral is taken over "differential 

spectrum" the resulting energy deposition values are smaller than A(z^) from Table 2a 

by factors ranging from 2.4 to 5.5. This is consistent with the energy (or 

ionization) -content discrepancy noted by Shemansky e_t al_ for an 

1.25 

incident electron spectrum with a T^ * power law. 


The distribution n^(z ) is not available from Maeda's work. We have therefore 
omitted the first term on the right-hand side of Eq(5) and have compensated for 
this by using as weighting factor the total stopping power L(T) rather than the 
restricted stopping power L^(T). The error incurred thereby is estimated to be 
small compared to the discrepancies by factors up to 5 which are being investigated. 

We have attempted, but failed, to account completely for this discrepancy, 
but were at least able to reduce it considerably, by analyzing what quantity was 
actually computed in Mae da's work.- The following points were considered: 

(1) The Monte Carlo program for simulating electron tracks used by Mae da 
was based on the procedures of BERGER (1963) and was thus an ancestor of th^ 

Carlo Model A used in the present work. The early version lacked certain refine- 
ments but. .was not essentially different from the current version. 

(2) The early version of the program included only primary electrons. How- 
ever, it can be seen from Fig. 8 that for the 20-keV source the contribution of 
secondary electrons becomes significant only at spectral energies ^ 4 keV. 

(3) The early version of the program provided as output not the electron 
flux spectrum in a semi-infinite medium, but the elec tron current transmitted 


**£n the Monte Carlo calculation the transmitted current is estimated by counting 
the electron tracks that cross a unit area of the exit surface of the slab tar- 
get. Once the electrons have left the slab they are not allowed to return to it, 
so that the exit surface is crossed at most once. The flux in a semi-infinite 
medium is estimating by computing the lengths of electron tracks in many thin 
layers of the medium parallel to the boundary plane. Each layer .of the medium 
is allowed to be crossed repeatedly by the same electron, in. the direction away 
from or toward the boundary plane. 


through a slab target. This current, denoted by Maeda as "differential energy 
spectrum" differs both conceptually and numerically from the electron flux 
spectrum as defined above. 
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In order to determine whether the discrepancy could be accounted for by 

the substitution of current for flux in Mae da's work, we have made model studies 

in which both quantities were 'calculated from the same set of sampled electron 

Monte Carlo histories. A typical set of results is shown in Fig. 10 which gives 

the total flux (denoted in this instance as F^T^)), the contribution to the 

flux from electrons moving in the direction of increasing depth, F 2■^ , ^* Z In^ , 

the current of electrons transmitted through a slab target of thickness z^, 

J(T z ). Of these three distributions, j(T,z ) is closest to Maeda's "differential 
x 1 m m 

energy spectrum", but is still higher by a factor^ of 2 or more. 

Some New Results . We would like to repeat that the spectra given here are 
preliminary, because they are calculated with oversimplified cross sections and by 
applying Monte Carlo Model A at lower energies than is, strictly speaking, per- 
missible. Furthermore it would be desirable to extend the spectra below the 
cut-off A of the present work of T q ) . ■ In any case, we have verified thau the 

energy content of the spectra is correct, in the sense that the relation between 
energy deposition and integrated spectra according to Eq(5) is satisfied (to 
within 5% at 20 keV, 5-10% at 10 and 5 keV, and 10-20% at 2 keV) .* 

ju 

*The decreasing accuracy at low source energies results from numerical approxi- 
mations of the Monte Carlo scheme. At 2 keV we have found it necessary to use 
the total rather than the restricted stopping power in Eq(5) to get the desired 
energy balance, corresponding to the fact that the Monte Carlo treatment at 
extremely low energies approaches the continuous-slowing— down approximation. 

The results to be shown indicate that the dependence of the electron flux 
spectra on the initial energy T q and the depth z -can be simplified considerably 
by certain scaling procedures. The dependence on can be minimized by expressing 
the spectra as functions of the ratio T/T , where. T is the spectral energy. The 
dependence on the depth can be reduced by plotting instead of F(T,z m > the 
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dimensionless quantity [t fv A(z ) 1 F(T z ) vs T/T , where r is the practical 
• op m o a o p 

range. The results in the following figures are normalized to one incident 
electron. 

The residual dependence of the scaled flux spectra on T q is illustrated in 

Fig. 'll. for T between 2 and 20 keV, and can be seen to be rather small. The 
,o 

flux spectra shown in the remaining figures all. pertain to electrons incident with 
energy T q = 10 keV. The change of the spectral shape with increasing depth is 
illustrated by Fig. 12. As expected, there is a spectral peak at or near T , due 
to primary electrons which have lost only a small part of their energy. This peak 
• . , disappears with increasing depth. There is also a low-energy component 

in the spectra, contributed (a) by primary electrons that have lost most of their 
energy in the course of traveling deep into the medium and were eventually turned 
around; (b) by low-energy secondary electrons resulting from inelastic collisions. 
Fig. 13 separates the flux spectra into the contributions from electrons traveling 
in the direction of increasing depth (DOWN) or decreasing depth (UP). At shallow 
depths one can distinguish three spectral regions. Near T q , downwar d-direc te d 
primary electrons predominate; at somewhat lower energies, upwar d-direc te d electrons 
predominate, and at low energies the number of upward and downward directed elec- 
trons is approximately equal, indicating that the low-energy flux component is more, 
or less isotropic. At great depths the downwar d-direc te d flux exceeds the 
the upward-directed flux, except at very low energies where the two components re- 
main approximately equal. Finally, Fig. 14 shows the dependence of the flux 
spectra on the direction of incidence, 0 q . This dependence is more pronounced 
for spectral energies near T^ than at low energies, because the low-energy elec- 
trons have undergone so many interac tions that they have "forgotten” their initial 
direction. 
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Table 1. Change of pitch angle, AS, due to the action of the geomagnetic 
field, for electrons with initial pitch angle 9^ traveling from 
height h^ to height h^. P (9 q ) is the percentage of electrons 
in a cosine-law distribution (IDH case) that have initial pitch- 
angles greater than 9 . 









<VV A(z m > 

Z m^ r p T , keV j Z n/ r 0 T , keV 

2 ____5_ 10_ 20 j p 20 50 100° 200 500 

0.0 j 0.68 0.62 0.61 0.61 j 

! I 

0.1 ; 0.97 0.89. 0.89 0.89 I 

0.2 j 1.17 1.10 1.12 1.13 I 


0.3 | 1.29 1.25 1.28 1.28 j 

i 

0.4 j 1.32 1.31 1.33 1.33 

i 

0.5 | 1.25 1.28 1.31 1.29 


0.6 j 1.09 1.13 1.16 1.15* | 

j 

0. 7 0.87 0.91 0.89 0.93 

0.8 0.58 0.60 0.59 0.63 

0.9 0.31 0.34 0.34 0.34 

i j 

! 1.0 0.15 0.18 0.16 0.16 I 


Table 2a. Scaled energy deposition function, for the case of a 
incident electron flux. 
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Table 2b. Scaled energy deposition function, for the case of an 
incident electron flux isotropic over the downward 


hemisphere 








T 

o 

(keV) 

r 

o • 

(g cm"' 2 ) 

a/ 

r 

P 

(g cm 2 ) 

b/ 

r 

P 

(g cm 2 ) 

2 

1. 86 (-5) 

1.52 (-5) 

1.40 (-5) 

5 

8.73(-5) 

7.19 (-5) 

6.65 (-5) 

10 

2.91(-4) 

2.48(-4) 

2 . 30 (-4) . 

20 

9. 82 (-4) 

8.45 (-4) 

7.83(-4) 

50 

4 .92 (-3) 



100 

1. 63(-2) 



200 

5.09 (-2) 



500 

2.00(-l) 



1000 

4.91(-1) 



2000 

i no" 




a/ For, .perpendicular incidence 
,b/ For the IDH case 

Table 3. Mean range, r^, and practical range, r^, for electrons in air. 



z 

m 

h 

Percent Contribution 

(g cm 2 ) 

(km) 

°<f = 5 keV 

e<=50 keV 

0 

300 

63 

29 

IQ’ 6 

194.6 

49 

19 

2 x 10 ^ 

! 

173.3 

40 

'15 

5 x 10 ^ 

1 

149.8 j 

21 

' 7 

10 -5 

135.7 ! 

7 j 

:l_I 


Table 4. Estimated contribution of electrons with initial energies T <^2 keV 

o 

to the energy deposition at various atmospheric depths. The 
incident electron flux is assumed to be isotropic over the 
downward hemisphere with an exponential energy spectrum 
proportional to exp (-T /a) . 




Table 5 . 


Height h at which the altitude profile of energy deposition peaks, 

and the corresponding value D(h)A(h) of the altitude profile, for 

the case of an incident electron flux isotropic over the downward 

hemisphere with a power-exponential energy spectrum proportional to 

o 

exp(-T^/a). Assumed geomagnetic latitude is - 90 . 


■#< 

height h, km 


a 

(keV) 

ay = 0 

2 

- 5 

10 

20 keV 

5 

108 

106 

104 

100 

96 

10 

100 

100 

98 

• 96 

93 

20 

93 

93 

92 

91 

90 

50 

85 

85 

85 

84 

83 

100 

80 

80 

80 

79 

79 

200 

73 

73 

73 

73 

73 


Altitude profile 

A . A 

of energy deposition, D(h)A(h), eV/cm 

a 

(keV) 

ay = 0 

2 

5 

10 

20 keV 

5 

1,29 (-3) 

2. 05 (-3) 

3. 37 (-3) 

6.07 (-3) 

1 . 22 (-2) 

10 

3.14 

3.98 

. 5.33 

7.82 

1.36 

20 

7.01 

7.93 

9.35 

1-19 (-2) 

1.74 

50 

1.91 (-2) 

2.00C-2) 

2.14(-2 ) 

2.38 

2.89 

100 

3.93 

4.02- 

4.15 

4.38 

4.84 

200 

7.56 

7.64 

7.77 

7.99 

8.43 



Table ^ . Dependence of the altitude profile of energy deposition, D(h)A(h), 

on the geomagnetic latitude X , for the case of an incident electron 

flux isotropic -over the downward hemisphere with a power-exponential 

energy spectrum proportional to T exp(-T /a), with y - 0.5 and 

^ per incident electron/cm 

a = 10 keV.. D(h)A(h) is given in units of eV/cm / Also shown is the 
height ft at which the altitude profile peaks. 


\ X 
\ m 

h T kriK 

o 

o 

cr» 

75° 

o 

o 

vD 

45° 

30° 

15° 

200 

1.73 (-5) 

1.73(-5) 

1.73 (-5) 

1.72(-5) 

1.70(-5) 

1.60(-5) 

180 

3.43 

3.43 ' 

3.42 

, 3.40 

3.35- 

3.09 

160 

7.79 

7.78 

7.75 

7.68 

7.50 

6.75 

150 

1 . 26 (-4 ) 

:..25(-4) 

1.25 (-4) 

1.24 (-4) 

1.20(-4) 

1.06(-4) 

140 

2.17 

2.17 

2.15 

2.13 

2.05 

1.77 

130 

4.20 

4.19 

4.16 

4.09 

3.91 

3.26 

i n A 
4. A.--VS 

9.93 

9.90 

,,9.80 

9.57 

9.00 

7.07 

110 

2.72 (-3) 

2.71 (-3) 

2.67 (-3) 

2.57 (-3 ) 

2. 34(-3) 

1.61 (-3) 

' 100 ■* 

5.23 

5.18 

5.03 

4.67 

3.87 

1.89 

90 

2.80 

2.75 

2.56 

2.17 

1.40 

2.89 (-4) 

85 

6.15 (-4) 

5.96 (-4) 

5.33(-4) 

4.10(-4) 

2.13 (-4) 

1.75 (-5) 

80 

3.33(-5) 

3.17 (-5) 

2.63 (-5) 

1.63 (* 5) 

4.61 (-6) 

- 

ft 

(km) 

98 

98 

98 

99 

100 

103 



O' 1 
-2 

g cm keV 

1*0 (-7) 
1 . 0 (- 6 ) 
2.0 
4 . 0 

1.0 (-5) 

2.0 

4.0 

1.0 (-4) 

2.0 

4.0 

l»0(-*3) 

2.0 

4.0 

1 . 0 (- 2 ) 

2.0 

4.0 

1.0 (- 1 ) 


20 

-1 

50 

100 

T , keV 
o 

.200 

500 

3.1(-4) 

5. 8(-4) 

7.0 (-4) 

5. 7 (-4) 

4.9(-4) 

3.1 

6.0 

7.0 

5.8 

5.0 

3.1 

6 . i 

7.0 

5.6 

5.0 

3.1 

6.2 

7.1 

5.8 

5.1 

3.2 

6.5 

7,2 

6.1 

5.2 

2.9 

6.8 

7.5 

6.5 

5.6 

2.5 

6.4 

7.9 

7.0 

6.2 

2.1 

5.0 

6.8 

7.1 

6.8 

1.8 

4.1 

5.1 

5.4 

6.2 

1.5 

3.2 

4.0 

3,8 

4.2 

1.1 

2.3 

2.6 

2.3 

2.4 

8. 3(-5) 

1.6 

1.7 

i.5 

1.5 

5.8 

1.0 

1.0 

• 8. 8 (-5) 

9.9 (-5) 

2.8 

4. 6(-5) 

4.5 (-5) 

4.3 

5.3 

1.2 

2.1 

2.3 

2.4 

2.6 

3. 8(-6) 

8. 5 (-6) 

1.1 

1.1 

7.6 (-6) 

4. 0(-7) 

2.6 (-8) 

7 *5 (-11) 

1.9 

3. 7 (-7) 
2.4 (-8) 

2.6 (-6) 
5. 4 (-7) 

2.0(-6> 



1000 

2000 

4.1(-4) 

4.1 (-4) 

4.2 

4.2 

4.3 

4.4 

4.4 

4.6 

4.8 

5.0 

5.2 

5.7 

5.8 

6.7 

6.6 

7.8 

6.7 

7.9 . 

4.9 

6, 3 

2.8 

4.0 

2.0 

2.9 

1.4 

1.9 

6. 4 (-5) 

5.8(-5) 


Table 7 . Energy deposition 
function for bremsstrahlu 
for the case of uniform 
wide-area precipitation 
of an electron flux 
isotropic over the down- 
ward hemisphere. The 

« 

quantity given is 

A, (z ) /T , in units 
dt m o 

2 -1 

of cm g . Numbers in 
parentheses indicate 
powers of ten. 
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Fig. 1 


Fig. 2 


Fig. 3 


Fig. 


Fig. 


Backs catering coefficients Rj, (number albedo) and Rj. (energy albedo) 
for monoenerge tic electrons incident on a semi-infinite air medium. 

The calculated coefficients are for perpendicular incidence^ = 0 ) 
or for an incident -flux isotropic over the downward hemisphere (XDH) . 

The experimental values of Rj, .are for perpendicular incidence and have 
been obtained by interpolating to atomic number Z - 7.4 the data 
the following authors: 

B: BISHOP (1966); H; HEINRICH (1966); VP: WEINRYb and 

PHILIBERT (1964) ♦ TV: TRUMP and VAN DE GRAAFF (1949); 

SA: SALDICK and ALLEN (1954); WT: WRIGHT and TRUMP (1962); 

GG: GLAZUNOV and GUGLYA (1964); S: STERNGLASS (1954); 

P; PALLUEL (1947); VA: VERDIER and ARNAL (1968). 

Backscattering of monoenerge tic electrons from a semi-infinite air 

medium, as a function of the incident obliquity angle 0 q . 

a. Number albedo, Ry 

b. Energy albedo, Rg 

Angular distribution W(9) of the current of backscattered 

electrons emerging from a semi-infinite air medium. The distributions . 

shown are normalized to unity over a 2rr-solid angle. 

^ Cumulative energy spectra of electrons backscattered from a semi-infinite 

air medium. The curves show, as a function of the incident energy T^, 

the fraction of the backscattered electrons that have energy T greater 

than 0.2 T , 0.4 T . 0.6 T or 0.8 T 
o o u y 

a. Perpendicular incidence, 9 q ” 0 

b. Incident flux isotropic over downward hemisphere, 0 q :IDH. 

5. Cumulative energy spectra of electrons backscattered from a semi-infinite 

air medium. The curves show, as a function of the incident obliquity 

angle 9 , the fraction of the backscattered electrons with energy T greater 

than 0.2T, 0.4 T , 0.6T or0.8T, where T is the incident energy. 

• q ’ q o o o 

a, T s 10 keV b. T q = 100 keV c. T Q = 1000 ^ eV 



Fig. 6. Altitude profile of energy deposition in the atmosphere, for an incident 
electron flux isotropic over the downward hemisphere with an exponential 
energy spectrum i exp(-T o /a). The relation between atmospheric depth and 
altitude is assumed to be that given by the CIRA (1965) Mean Atmosphere 
at geomagnetic latitude = 90°. Results are normalized to one incident 
electron/ cm 2 . 

Fig, 7. Altitude profile of energy deposition by electron-produced bremsstrahlung 
in the atmosphere, for the case of wide-area precipitation of an electron 
flux incident isotropically over the downward hemisphere with an exponential 
spectrum - exp (-T Q /a) . The relation between atmospheric depth and alti- 
tude is assumed to be that given by the CIRA (1965) Mean Atmosphere at 
geomagnetic latitude = 90°. Results are normalized to one incident 
electron/cm 2 . Dashed portions of curves are extrapolations. 


Fig. 8. Comparison of the energy flux spectrum F ( T » z m ) with Mae da's "differential 
energy distribution", for 20-keV electrons incident perpendicularly on an 
air medivm. The dotted curves indicate the contribution to F(T,z^} from 
degraded primary electrons. 

Fig. 9. Depth distribution of electrons at cut-off energy A. The ordinate is 
the dimensionless quantity n^(z^) , where ^- s number of 

electrons per unit depth whose energy falls below A, P and S indicate 
primary and secondary electrons, respectively. The normalization cor- 
responds to one incident primary electron. 

Fig .10. Comparison of various spectral distributions, at a depth z^ = 0.5 r Q> for 
20-keV electrons incident perpendicularly on an air medium. 

F, (T.z ): total flux spectrum 
4rr m . 

F (T z ) : contribution to total flux spectrum from electrons 
2rr v ’ m 


moving in the direction of increasing depth 

J (T.z ) : spectrum of the current of electrons transmitted 

m 

through a slab target with a thickness equal to z . 

Fig. 11. Dependence of the electron flux spectrum F(T,z^) on the incident energy T q . 



Fig- 12- 


Dependence of the electron flux spectrum 
incident with an energy of 10 keV. 
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Fig. 13- Directional characteristics of the electron flux spectrum, for electrons 
incident with an energy of 10 keV. The curves labelled DOWN and UP 
. , represent.; the contrive tivhs ...from ..elec Lrous ..moving In the direction of 
increasing and decreasing depth, ...respectively. The -curve marked TOTAL is 

“the sum of DOWN and UP. 

Fig. l4. Dependence of the flux spectrum F(T,z ) on the obliquity (pitch angle) 6 
of the incident electrons, for an incident electron energy of 10 keV; 
the scaling parameter in this figure is the practical range for the 
case of perpendicular incidence. 
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